Spectral Analysis of Excited Nucleons in 
Lattice QCD with Maximum Entropy Method 



Kiyoshi Sasaki a , Shoichi Sasaki b,c , Tetsuo Hatsuda b , 

a Center for Computational Science, University of Tsukuba, Tsukuba, Ibaraki 

305-8577, Japan 

h 'Department of Physics, University of Tokyo, Tokyo 113-0033, Japan 

C RIKEN BNL Research Center, Brookhaven National Laboratory, Upton, NY 

11973-5000, USA 



Abstract 

We study the mass spectra of excited baryons with the use of the lattice QCD 
simulations. We focus our attention on the problem of the level ordering between the 
positive-parity excited state iV'(i440) (the Roper resonance) and the negative-parity 
excited state iV*(1535). Nearly perfect parity projection is accomplished by combin- 
ing the quark propagators with periodic and anti-periodic boundary conditions in 
the temporal direction. Then we extract the spectral functions from the lattice data 
by utilizing the maximum entropy method. We observe that the masses of the N' 
and N* states are close for wide range of the quark masses (M n = 0.61 — 1.22 GeV) 
and for the physical point after the extrapolation in terms of the quark mass. The 
latter aspect is in contrast to the phenomenological prediction of the phenomenolog- 
ical quark models. The role of the Wilson doublers in the baryonic spectral functions 
is also studied. 



1 Introduction 



The lattice QCD simulations have been shown to be a powerful tool to study 
the properties of hadrons from first principle. In particular, ground state 
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hadron masses in the quenched approximation agree with the experimental 
values in 10% accuracy [1]. Furthermore unquenching of the dynamical quarks 
tends to improve the agreement [2]. On the other hand, studies of the excited 
hadrons are still in an exploratory stage in spite of their physical importance. 
(For a brief review, see [3] and references therein.) 

In this paper, from the quenched lattice QCD simulations with Wilson fermions, 
we study excited nucleons with a special focus on the spectrum in the positive- 
parity channel A/ 7 (1440) (the Roper resonance) and that in the negative-parity 
channel iV*(1535). 

The fact that the mass of the N'(14A0) is smaller than that of the A/"* (1535) has 
been a long standing puzzle: Indeed, phenomenological quark models predict 
that M N , - M N ~ 2 x (M N * - M N ) because N = (Is) 3 while N' = (ls) 2 (lp) 
and N* = (ls) 2 (2s) or (ls)(lp) 2 [4]. As possible resolutions of this puzzle, 
exotic descriptions of the N' state and/or the iV* state have been proposed 
such as the chiral quartet scheme [5] and the pentaquark picture [6] as well as 
the quark model with the effective flavor-spin interaction [7]. 

In principle, the lattice QCD simulation ought to give the final answer to this 
problem. However, definite conclusions have not been drawn yet even in the 
quenched lattice QCD studies because of several technical reasons: First of all, 
one needs to make appropriate projection and isolation of the excited states 
to differentiate the N, N' and N* states. Furthermore, lattice artifact due to 
finite lattice volume becomes severer for the excited states than the ground 
state [8,9]. 

To overcome these problems, (i) we adopt a technique to combine the quark 
propagators with periodic and anti-periodic boundary conditions in time for 
precise parity projection [10], (ii) we employ the maximum entropy method 
(MEM) [11,8] for extracting the excited mass spectra, and (iii) we utilize 
a larger spatial size (La ~ 3.0 fm) with the lattice spacing a ~ 0.093 fm 
03 = 6.0). 

Our main conclusion is that N' and iV* states are nearly degenerate within 
statistical errors in the wide range of the quark mass (M w = 0.61 — 1.22 GeV). 
We found that, after the chiral extrapolation with an empirical "curve fit" 
formula, the masses of the N' and iV* states are consistent with the approx- 
imate degeneracy observed in experiment. Also, we found a firm evidence of 
high-lying bound states composed of Wilson doublers. 

This paper is organized as follows. In Section 2, we provide a brief introduc- 
tion to the parity projection and MEM. After summarizing our simulation 
parameters in Section 3, we discuss the numerical results obtained from the 
MEM analysis of the spectral function in Section 4. The bound states of Wil- 
son doublers which appear as high energy peaks in the spectral function are 
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discussed in Section 5. A brief comparison with previous studies is given in 
Section 6. Section 7 is devoted to summary and concluding remarks. 



2 Parity projection and MEM 



The nucleon is described by the trilinear quark composite operator given by 



N(x) = e abc [u^(x)C-f 5 d b (x)]u c (x), (1) 

which is composed of a local scalar-diquark operator and a spectator-like quark 
field. This local composite operator couples to both positive and negative 
parity states [12,10]. The asymptotic form of the zero-momentum nucleon 
correlation, G(t), may be expressed in finite temporal extent, T, as 



G(t) = ]T(A(f,t)A(C),0)> 

X 

— > [k + A + e~ M+t - A_A_e" M -*" 

+b [A_A + e~ M +^ - A + A_e- M -( T -V] (2) 

with A± = (1 ± 74) /2. The coefficient b = 1(— 1) is for the periodic (anti- 
periodic) boundary condition in the temporal direction. M + (M_) is the ground 
state mass of the positive (negative) parity nucleon. The second and forth 
terms on the right hand side of Eq.(2) appear since the nucleon is a composite 
particle. Due to the presence of the third and forth terms corresponding to the 
primal reflections from the time boundary, the naive projection Tr (A±G(t)) is 
insufficient to extract a contribution of the desired parity state. 1 A solution 
to this problem is to realize the free-boundary situation (6 = 0), which can be 
achieved by using a linear combination of the quark propagators with periodic 
and anti-periodic boundary conditions in time [9,10]. The linear combination 
in the quark level automatically fulfills the linear combination in the hadronic 
level and then manages the free-boundary situation up to the first wrap-round 
effect [9]. Adopting this procedure, one can project out the desired parity state 
by A± accurately without the contamination from the opposite parity state. 

The baryonic correlation G(t) with the above parity projection has a spectral 
representation, 

G(t) = / duo K(t,u)A(u), (3) 
Jo 



1 The most primitive prescription to avoid the contamination from the opposite 
parity states is to analyze the data only up to about T/2. 



3 



where K(t,u>) = exp(-ujt) is the Laplace integral kernel. The spectral func- 
tion (SPF), A(u>), carries all the information of the hadronic states (both 
ground and excited states). The discretized version of Eq. (3) reads G(tA = 
J2iexp(—LUiti)A(uJi): We need values of A{uo{) for sufficient number of frequen- 
cies ui, typically of (9(100). The Monte Carlo data G(U) is, however, available 
only at some dozens of time slices on the lattice. This is a typical ill-posed 
problem, where the number of equations is much smaller than the number 
of unknown quantities. The maximum entropy method (MEM) introduced in 
Ref. [11] enables us to determine SPF uniquely by supplementing the lattice 
data with prior knowledge such as the positivity and the asymptotic behavior 
of SPF. 

The conditional probability of A(u) given the Monte Carlo data G(t) and the 
prior knowledge H is given by P[A\GH] which is evaluated by using the Bayes' 
theorem as 

P[A\GH] oc P[G\AH] P[A\H}. (4) 

Here P[Cr|Aff] is the standard likelihood function and PL4|if] is given by the 
Shannon- Jaynes entropy: 



P[A\H] oc exp 



1=1 L 



Given G(t), a and A, the most probably A is obtained by the stationary 
condition, (8 /8A)P[A\GH\ = 0, which has a unique solution if it exists [11]. 
Also, the reliability of the solution A can be tested by the second variation, 
(5/5A) 2 P[A\GH}. We have A = m cu n (n = 2(5) for meson (baryon)) in 
Eq.(5): This is called the default model in which the parameter mo may be 
fixed either by the perturbative asymptotic behavior of SPF at large uj, or 
by demanding that the resultant A has smallest error bar. By evaluating the 
probability P[a\GH] using the Bayes' theorem, final form of the spectral func- 
tion is obtained by averaging over a with the weight factor P[a\GH]. For more 
details on MEM techniques, see Ref. [11]. 



3 Simulation parameters 



We have performed simulations in quenched QCD with the single plaquette 
gauge action and the Wilson fermion action. We utilized the Metropolis algo- 
rithm with 20 hits at each link update. The gauge ensembles in each simulation 
are separated by 1000 sweeps for (3 = 6.0, and 600 sweeps for (3 = 5.8. For 
the matrix inversion, we used BiCGStab algorithm [13] with the convergence 
condition |r| < 10~ 8 for residues. We calculate the point-to-point nucleon 
propagator constructed by the point-to-point quark propagator so that the 
spectral function is assured to be positive semi-definite. Our calculation has 
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Table 1 

Simulation parameters in this study. (3 = 6/g 2 , L (T) is the number of the spatial 
(temporal) site, a (La) is the lattice spacing (spatial lattice size) in the physical unit, 
and iV con f is the number of gauge configurations. The Sommer parameter r$ = 0.5 
fm is used to fix the scale [15]. 



p 


L 3 x T 


a (fm) 


La (fm) 


hopping parameter 


-^conf 


5.8 


24 3 x 32 


0.136 


3.3 


0.1600, 0.1570, 0.1555, 0.1540, 0.1530 


200 


6.0 


16 3 x 32 


0.093 


1.5 


0.1550, 0.1530, 0.1515 


580 




32 3 x 32 


0.093 


3.0 


0.1550, 0.1530, 0.1515, 0.1500 


200 



been performed on a Hitachi SR8000 parallel computer at KEK (High Energy 
Accelerator Research Organization), using the extended code based on the 
Lattice Tool Kit (LTK)[14]. 

Our simulation parameters are listed in Table 1. Physical scale is set by the 
Sommer parameter r = 0.5 fm [15]. Our main conclusion will be drawn from 
the data taken on the finest and largest lattice (f3 = 6.0 and L = 32) in the 
Table 2. The case with coarser lattice spacing (f3 = 5.8 and L = 24) is also 
studied to analyze the nature of the high-lying peaks of the spectral function 
associated with the Wilson doublers. A supplementary data with a smaller 
lattice size (L = 16 with (3 = 6.0) are also taken to check the effect of finite 
lattice size on the excited nucleon spectrum. For each data, several hopping 
parameters are used to cover the range M n ~ 0.6 — 1.2 GeV. 



4 iV*(1535) and the Roper resonance iV'(1440) 

In Fig.l, we have shown examples of the correlation functions G(t) in the 
positive parity and negative channels for (3 = 6.0 (a = 0.093 fm), k = 0.1550 
(M n ~ 0.6 GeV) and L = 32 (La ~ 3.0 fm). Fig.2 show the dimensionless 
spectral functions, p(uo) = A(uj)/uj 5 , obtained from the MEM analysis of G(t) 
in Fig.l. The default model in the MEM analysis for the nucleon channel is 
chosen to be about mo = 0.004, which is deduced from the asymptotic form of 
the SPF at the one- loop level in the perturbative QCD in the continuum [8]. 
We have check that the final results are insensitive to the choice of m around 
this value. The temporal interval of the Monte Carlo data G(t), used in the 
MEM analysis, is chosen to be [t min ,t max ] = [1,20]. Exceptions are for the 
negative-parity nucleon with light quark masses (k = 0.1530, 0.1550) in which 
^max — 19 is used because of the noisy signals for t > 20. Note that t — 
corresponds to a source location and is excluded in the analysis. It is worth 
mentioning that sufficiently large statistics makes the reconstructed image of 
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Fig. 1. Examples of the correlation function G(t) generated by the Monte Carlo 
simulations as a function of Euclidean time t for f3 = 6.0, k = 0.1550 and L = 32. 
The left panel corresponds to the positive-parity nucleon channel and the right panel 
to the negative-parity nucleon channel. In the right panel, we have not shown the 
data for t > 20, since the statistic errors become much larger than the mean values 
and the data do not have statistical significance. 



SPF stable against variation of t max beyond some critical point around t max ~ 
18 — 19. In our MEM analysis, we have chosen the region, 0.2/ a < uj < 2n/a, to 
cover all the relevant states available on our lattice. We have carefully checked 
that changing the lower (upper) limit of u to smaller (larger) values have no 
effect on the finial results. 



The crossbars attached to the SPFs in Fig. 2 are attached with the procedure 
given in [11]. First of all, we arbitrarily choose an interval in the frequency 
space so that it covers one of the peaks in SPF. The horizontal bar in the figure 
represent this interval and is not at all related to the uncertainly of the peak 
position. After setting the interval, the uncertainty of the peak height averaged 
over the interval can be calculated from the second variation of P[A|Gif] 
with respect to A. The vertical bar in the figure represents ±ler error of such 
averaged peak height. Fig. 2 shows that the peaks are statistically significant 
except for the marginal second peak in the lower panel. We identify the first 
(second) peak of the upper panel in the figure as N (N'), and the first peak 
in the lower panel as N*. 

Carrying out the MEM analysis in the meson and baryon channels for different 
values of k and identifying the positions of low-lying peaks, we extracted the 
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Fig. 2. The dimensionless spectral function, p(uj) = A(lj)/uj 5 , as a function of the 
dimensionless frequency wa for (3 = 6.0, k = 0.1550 and L = 32. The upper panel 
corresponds to the positive-parity nucleon channel and the lower panel to the neg- 
ative-parity nucleon channel. 



hadron masses as summarized in Table 2. The statistical errors in the paren- 
theses are evaluated by the jack-knife method. Critical hopping parameter k c 
at which the pion mass vanishes reads k c = 0.15668(21) for (3 = 6.0. Our sim- 
ulations are done for quark masses which are not very close to the chiral limit 
so that a guidance by the chiral perturbation theory is no longer useful for the 
extrapolation of hadron masses down to the physical point. Furthermore, the 
leading chiral behavior of the excited baryon masses is not known at present. 
In this situation, adopting a simple formula, which can describe the observed 
dependence of the quark mass well, is a reasonable choice to make the chiral 
extrapolation. In this paper, we utilize the "curve fit" ; 

(aM H ) 2 = do + d 2 (aM w ) 2 . (6) 

which does not include a term linear in aM n , which is responsible for the 
leading chiral behavior for the nucleon in quenched chiral perturbation the- 
ory [16]. As reported in Refs. [17,9], Eq.(6) gives a better fit to the data in 
the heavy-quark region than the linear fit. Note that we have not made ex- 
haustive study on the chiral extrapolation in this paper, since the systematic 
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errors stemming from the finite-size effect is more significant than those from 
the chiral extrapolation. 

As for the ground state in each quantum number, one may compare the MEM 
result and the conventional single exponential fit using G(t — > large). Shown 
in Fig. 3 are the comparisons of the masses of the iV (the ground state in 
the positive-parity nucleon) and the N* (ground state in the negative-parity 
nucleon) in two approaches: The circles represent the mass evaluated by the 
single exponential fit, while the squares represent the position of the first peak 
obtained by the MEM analysis. They are consistent with each other within 
the statistical error for the whole range of the quark masses. 

Let us consider the finite size effect on the peaks of the spectral functions. In 
Fig. 4, we compare the SPFs on a larger lattice (La ~ 3.0 fm) shown by the 
solid lines with those on a smaller lattice (La ~ 1.5 fm) shown by the dashed 
lines. Although the position of the N (the first peak in the upper panel) has a 
minor shift, the positions of the N' (the second peak in the upper panel) and 
the ./V* (the first peak in the lower panel) have significant upward shifts as we 
decrease the lattice volume. This is consistent with an intuitive observation 
that the size of excited hadron is larger than that of the ground state and thus 
receives a large finite volume effect. Detailed analysis on the finite size effect 
on the N and the iV* by two of the present authors shows that the spatial size 
of about 3 fm is necessary to obtain reliable results even for relatively heavy 
quark masses (M w ~ 1 GeV) [9]. With these observations, we will use the data 
for La ~ 3.0 fm in this paper for final results. 

Now, we compare the position of the Roper resonance (N') and the negative 
parity nucleon (N*). Plotted in Fig. 5 are the masses of these resonances as well 
as the nucleon mass as a function of the pion mass squared in the physical unit. 

Table 2 

Masses of the pion, the p-meson, the nucleon N (J p = l/2 + ), the Roper N' (J p = 
l/2 + ) and the N* (J p = l/2~) extracted from the data on the lattice with L = 32 
and (5 = 6.0 by the MEM analysis. The statistical errors are evaluated by the 
jack-knife method. The values in the chiral limit are evaluated by the curve fit. 



K 


aM n 


aM p 


aM N 






K c 





0.318(15) 


0.464(21) 


0.76(16) 


0.689(102) 


0.1550 


0.288(14) 


0.422(13) 


0.635(14) 


0.92(13) 


0.843(68) 


0.1530 


0.420(7) 


0.505(4) 


0.791(9) 


1.09(8) 


1.030(56) 


0.1515 


0.502(3) 


0.568(3) 


0.898(7) 


1.20(6) 


1.143(50) 


0.1500 


0.578(2) 


0.630(5) 


1.000(6) 


1.31(6) 


1.236(49) 
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Fig. 3. A comparison between the position of the first peak in the SPF obtained by 
the MEM analysis and the mass evaluated by the single exponential fit, for L = 32 
and P = 6.0. The left panel correspond to the positive-parity nucleon channel and 
the right panel to the negative parity nucleon channel. 



All the points are evaluated by the MEM analysis. The chiral extrapolations 
with the curve fit are also shown together with the experimental values. We 
find that the N' and N* states are almost degenerate within statistical errors 
in a wide range of quark masses, 0.61 < M.,, < 1.22 GeV, although the central 
values of the N' is slightly higher than the N* in this interval. In the chiral 
limit, we have the following mass ratios: 

M N ,/M N = 1.64(35) (Expt. ~ 1.53), 

M N */M N = 1.48(22) (Expt. ~ 1.63). (7) 

Our result here is in contrast to what is expected in the phenomenological 
quark models as mentioned in the Introduction, but is consistent with the 
experimental value. To make a firm conclusion on the precise level ordering 
between the N' and N* states, we need to increase the number of gauge con- 
figurations to reduce the statistical errors and also need to make simulations 
with smaller quark masses. 

In relation to the physics in the vicinity of the chiral limit, we comment on 
the quenched artifact due to the rj'N ghost state in the nucleon correlator. It 
is pointed out in [18] that this ghost state affects the nucleon correlator sig- 
nificantly for ^ 0.3 GeV. If such effect is significant, the spectral function 
becomes negative and the standard MEM analysis assuming the positivity of 
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Fig. 4. The finite size effect on the SPF for (3 = 6.0 and n = 0.1550. The upper 
(lower) panel corresponds to the positive (negative) parity nucleon channel. In each 
panel, we compare the SPF obtained from the larger size, L = 2>2{La ~ 3.0 fm), 
with that from the smaller size, L = 16(La ~ 1.5 fm). Note that the 2nd lowest 
peak on L = 32 lattice disappears on L = 16 lattice in the negative parity channel. 



the SPF is invalidated. In the present study, however, simulations are per- 
formed in the region M n = 0.6 — 1.2 GeV, and we have checked explicitly that 
the rj'N threshold is higher than the masses of the N, N' and N* states in the 
this interval. 



5 Bound state of doublers 



In this section, we discuss the high-lying peaks of the SPF which could be the 
bound states containing Wilson doublers. This possibility was first pointed out 
by CP-PACS collaboration for the meson spectra [19]. Since the mass of the 
Wilson doublers is inversely proportional to the lattice spacing a, studying the 
peak positions of SPF as a function of ua provides us with a useful information 
on the doubler bound states. 

Shown in the left panel of Fig. 6 are the SPFs in the positive parity nucleon 
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Fig. 5. Ground and excited nucleon spectra as a function of the pion mass squared 
in the physical unit. The experimental values for iV(940), iV'(1440) and iV*(1535) 
are also marked by filled diamonds. 



channel for M w ~ 1 GeV with (5 = 5.8 (upper figure) and f3 = 6.0 (lower 
figure). The lattice volume is taken to be La ~ 3.0 fm. The MEM analysis is 
performed with t max = 14 for (3 = 5.8 and 20 for (3 = 6.0. These values of t max 
are chosen so that the resulting images of SPF are stable enough. The number 
of configurations in both cases is iV con f = 200. 

Consider the two high-lying peaks whose positions in terms of ua are almost 
independent of the lattice spacing a. Insensitivity of the peak positions in 
terms of ua is already an indication that they contain Wilson doublers, the 
pole mass of a free Wilson doubler is given by am^ c = \n(l + am q + 2n) with n 
being an integer. The bare quark-mass m q is defined by am q = — l//c c )/2. 
Note that baryons at rest composed of Wilson doublers need to contain two 
doublers to satisfy the momentum conservation. Then one may write a simple 
mass formula for two high-lying peaks as 



aM„^o = 2 ln(l + am q + 2n) + m(l + am q ) + C n , (8) 

where C n corresponds to a binding energy between the two doublers and a 
physical quark. As far as n ^ 0, aM n depends only weakly on the bare quark- 
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Fig. 6. Left panel: A comparison of the SPF in the positive parity channel at two 
different couplings, [3 = 5.8, 6.0. The hopping parameter is taken as k = 0.1530 
(k = 0.1515) for = 5.8 (f3 = 6.0) so that we have M w ~ 1 GeV in both cases. Right 
panel: The positions of the two high-lying peaks, labeled by 1st UPS (first unphysical 
state) and 2nd UPS, are plotted against the bare quark-mass am q = (1/re — 1/k c )/2. 
Results of fitting the data using a mass formula given in the text are shown by the 
dotted lines. 



mass m q . 

In the right panel of Fig. 6, the positions of the highest two peaks for j3 = 6.0 
are plotted as a function of the bare quark- mass, am q . Results of fitting the 
data by Eq. (8) are shown by the dotted lines. We find that the mass formula 
reproduces the positions and the m q dependence of the highest two peaks 
well. As a result of this quantitative analysis, it is quite likely that the highest 
two peaks are the bound states containing Wilson doublers which eventually 
decouple from the physical spectrum in the continuum limit, a — > 0. 

It is worth mentioning that any prominent structure of SPF in the higher 
frequency region {u ^> 1/a) is not found in the 4-dimensional 0(4) a model 
[20] and the 3-dimensional four-fermion model with staggered fermions [21]. 
These observations support that our observed high-lying peaks far beyond 
cutoff scale 1/a are associated with Wilson doublers. 



6 Comparison to other approaches 



In this section, we briefly make a comparison among different approaches ap- 
plied to analyze the excited nucleons. (See Table 3.) The methods to extract 
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Table 3 

Summary of the simulation parameters and methods in previous analyses together 
with the current analysis. The fermion action, the lattice spacing, the pion mass, the 
spatial size, the boundary condition in time and the method of analysis are tabu- 
lated. CID stands for the chirally improved Dirac operator, CCF for the constrained 
curve fitting. Note that the physical scales are set by the pion-decay constant f w in 
Ref. [18], and by the Sommer parameter vq in other references. 



fermion 


a- 1 (GeV) 


M w 


(GeV) 


La (fm) 


B.C. 


method (Ref.) 


DWF 


2.1 


0.57 


- 1.43 


1.5 


P+AP 


2x2 correlator [10] 


clover 


2.9 


0.51 


- 1.08 


1.6 


P 


modified correlator [22] 


CID 


1.3 


0.45 


- 0.87 


1.8 


AP 


4x4 correlator [23] 


overlap 


1.0 


0.18 


-0.87 


2.4, 3.2 


AP 


CCF [18] 


Wilson 


2.1 


0.61 


- 1.22 


1.5, 3.0 


P+AP 


MEM (this work) 



the excited spectra may be classified into three categories: (i) Diagonalization 
method using correlation matrix [10,23], (ii) Subtraction method using modi- 
fied correlator [22] and (iii) Bayesian methods such as the maximum entropy 
method ([8] and the present work) and the constrained curve fitting (CCF) 
[18]. 

In (i), one starts with the correlation matrix, CV,(t) = (Oj(t)Oj(O)), and solve 
the generalized eigenvalue problem C(t)v = X(t, t )C(t )v where t is fixed and 
t > t . The k-th largest eigenvalue of the transfer matrix X(t, to) may be given 
by A^(t, t ) = e~('~'°) Mfe with appropriate choices of t and t [24]. In Ref. [10], 
an unconventional operator, e a bc[ul{x)Cdb{x)]^u c {x) is utilized to construct 
the 2x2 correlation matrix in addition to the conventional one as Eq.(l). On 
the other hand, in Ref. [23], they build up the correlation matrix with four 
operators, which correspond to four possible combinations to construct the 
nucleon operator from two sets of the smeared quark propagator. In (ii), a 
modified correlator, G{t) = G(t+ l)G(t — 1) — G 2 (t), is investigated instead of 
the original correlator G(t), so that the sum of the first two states, M + M 1: 
governs the leading-exponential decay of the modified correlator as G(t) = 
2E„> A)^n[cosh(M„ - M ) - 1] e -( Mo+M ")* [22]. In (iii), prior knowledge of 
the correlator is considered in addition to the Monte Carlo data. In MEM 
[11], the prior knowledge is introduced in the entropic form, while in CCF, it 
is introduced in the gaussian form [25]. 

In Refs. [10,23,22] (the first three works in Table 3), degeneracy between the 
N' and N* states has not been observed within the range of the quark masses 
where simulations are performed. On the other hand, in Ref. [18] and the 
present work (the last two works in Table 3), an approximate degeneracy 
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of the N' and the iV* are seen within the statistical errors, although it is not 
yet clear where and how the level switching between the N' and N* states 
really takes place as a function of the quark mass in our study. One of the 
major differences between the first three works and the last two works is the 
spatial lattice size. 

As we have discussed in Section 4 and will be reported in detail in a separate 
publication [26], our simulation with a smaller lattice (La ~ 1.5fm) show that 
there is a sizable upward shift especially in the N' channel for M n ~ 0.6 and 
0.9 GeV. This lattice artifact increases the N' — iV* difference even for rather 
heavy quarks and may partly explain the reason why degeneracy have not 
seen in first three works. 



7 Conclusion 

In this article, we have studied the mass spectra of excited nucleons in the 
quenched lattice QCD simulations. We have focused our attention on the level 
ordering between A/ 7 (1440) (the nucleon excited state in the positive-parity 
channel) and AT* (1535) (the nucleon excited state in the negative-parity chan- 
nel) which has been a long standing puzzle. Indeed, it is difficult to be ex- 
plained by phenomenological quark models in which the Roper resonance N' 
has much higher mass than the iV*. 

To attach the problem with a firm footing, we have utilized a technique to 
combine quark propagators with periodic and anti-periodic boundary condi- 
tions in time to eliminate the contamination of the opposite parity states in 
the baryonic correlator up to the first wrap-round effect. Also, we have used 
the maximum entropy method (MEM) to reconstruct the spectral function 
which contains not only the information of the ground state but also that of 
the excited states. Since the finite size effect is shown to be severe for excited 
nucleons [9], we have taken a large lattice size (La ~ 3 fm with (5 = 6.0). 

With all the above, we have extracted the masses of the N, N' and N* states 
for M n = 0.61 — 1.22 GeV, which are summarized in Fig.5. Our results show 
that the masses of the N' and A^* states are approximately degenerate within 
statistical errors for wide range of quark masses, although the central value of 
the N' is slightly higher than the A^*. This tendency is valid also at the physical 
point after making chiral extrapolation with a curve fit, Eq.(7). The latter 
result is in contrast to the predictions of the phenomenological quark models 
but is consistent with similar finding in Ref. [18] which employs different lattice 
action and different way of data analysis. Whether level switching between the 
N' and N* states takes place at certain quark mass is an open question left 
for future studies. We have also studied the high-lying peaks of the SPF by 
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carrying out simulations on a coarse lattice {(3 = 5.8) with approximately the 
same lattice volume. From the analysis of the a-dependence and the quark- 
mass dependence of these peaks, we found a firm evidence that the high-lying 
peaks correspond to the bound states composed of two Wilson doublers and 
one normal quark. 

There are various directions to be explored in future studies: To increase the 
statistics to pin down the precise location of the level switching; To generalize 
the MEM to treat the negative SPF originating from the ghost rj'N state; 
To make systematic studies of the other excited baryons and mesons. The 
nucleon excited states at finite temperature are also interesting to be studied 
in connection with realization of chiral symmetry in baryons [27]. The MEM 
analysis also may be helpful to disentangle the pentaquark signal from a tower 
of the KN scattering states [28]. 
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